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ABSTRACT 

We use Hubble Space Telescope archival images to measure central surface brightness profiles of globular 
clusters around satellite galaxies of the Milky Way. We report results for 21 clusters around the LMC, 5 around 
the SMC, and 4 around the Fornax dwarf galaxy. The profiles are obtained using a recently developed technique 
based on measuring integrated light, which is tested on an extensive simulated dataset. Our results show that 
for 70% of the sample, the central photometric points of our profiles are brighter than previous measurements 
using star counts with deviations as large as 2 mag/acrsec^. About 40% of the objects have central profiles 
deviating from a flat central core, with central logarithmic slopes continuously distributed between -0.2 and 
-1.2. These results are compared with those found for a sample of Galactic clusters using the same method. 
We confirm the known correlation in which younger clusters tend to have smaller core radii, and we find that 
they also have brighter central surface brightness values. This seems to indicate that globular clusters might be 
born relatively concentrated, and that a profile with extended flat cores might not be the ideal choice for initial 
profiles in theoretical models. 

Subject headings: globular clusters: general, stellar dynamics 



1. EVITRODUCTION 

The observational study of internal dynamics of glob- 
ular clusters (GCs) has benefited from imaging from 
space as well as enhanced sp ectroscopic capabilities in the 
ground. Novola & Gebharda ('2006) (From now on called 
"Paper I") measure surface brightness profiles (SB) from 
Hubble Space Telescope (HST) images for a sample of 38 
galactic globular clusters. The results from that work show 
that half of the objects in the sample are not consistent with 
having central flat cores, but instead, the distribution central 
surface brightness logarithmic slopes is continuous form -0.2 
to -0.8. The ages of the Galactic cl usters are all confined 
to a narrow range older than '-^lOGyr (' Salaris & Weissll2002l: 



'angi 

iDe Angeh et al.ll2005l) . It is desirable to measure central SB 
profiles of globular clusters with younger populations to find 
out if these central cusps are also observed in less evolved 
clusters. Globular clusters around Milky Way satellites are 
ideal targets for this task since they have a larger age range, 
they are relatively near, and many of them have been observed 
with HST. 

Surface brightness profiles have been obtained for GCs 
in the Large Magellanic Cloud (LMC), Small Magellanic 
Cloud (SMC), and Fornax dwarf galaxies in various stud- 
ies usin g ground-based data . For the LMC cluster s, star 
counts (Kontizas et al.l IT987h . aperture photometry (Mateo' 
[1987; Elson 1991) and hybrid techniques (Elsonet al. 1987) 
have been used to obtain surface density profiles for a va- 
riety of subgroups (rich, old, young, disk, and halo clus- 
ters). For the SMC clusters, only a fe w studies have mea- 
sured density profiles from star counts dKontizas & KontizasI 
[1983: Kontizas et al. 1986). A couple of studies measure 
density profiles from aperture photometry for globular clus- 
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ters around the Fornax dwarf galaxy dSmith et al.l 119961: 
!Rodgers & Roberts 1994). All of these studies are very use- 
ful for studying SB profiles at large radius, but at small radius 
they suffer from the usual seeing and crowding problems as- 
sociated with ground-based observations. 

A large systematic study of surface brightness profiles ob- 
tained from space-based imaging has been carried out by 
Mackey & Gilmore (2003a,b,c) (from now on collectively 
referred to as MAC03). They gather a broad sample of 
LMC, SMC and Fornax galaxy GCs imaged with WFPC2. 
They obtain SB profiles by measuring star counts weighted 
by brightness from which they derive fundamental quanti- 
ties l ike central density and core rad ius by fitting EFF pro- 
files (lElson et all lm% Ulsonl [T99l1) . which are power-law 
plus core profiles with three parameters: core radius, cen- 
tral surface brightness, and slope of the power-law. They de- 
termine that 20±7% of the clusters in their sample are con- 
sistent with a post-core-collapse morphol ogy, a simila r num- 
ber to the one found for Galactic clusters (iTrager et al.l ll995). 
When they compare their profiles with previous results ob- 
tained from ground based images, they find that important 
aspects of the nature of the profiles can be measured by im- 
proving the spatial resolution. McLaughlin & van der Mare| 
(12005) (from now on called MVM05) combine the data from 
MAC03 with star counts profiles from ground-based data in 
order to obtain a more accurate photometric normalization. 
They fit the re-normalized de-reddened resulting profiles with 
variety of models such as KingJits^(King 1966), an alternate 
modified isothermal model by lWilsonI d 19751) . which has more 
extended envelopes than a King model, and a power-law plus 
core model like the one used in MAC03. They conclude that 
The Wilson fits provide the best description of the outer part 
of the clusters for both o ld and young populations. 

son et all (Il989h and lElsonI ill 9921) find an interesting rela- 
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TABLE 1 
Simulation input. 
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input slope 

Fig. 1 . — Compaiison between the measured and input center for every set 
of simulations with different central slopes. The average distance between 
the actual and measured center is shown for the x (solid points) and y (open 
points) coordinates in pixels. A small horizontal offset is introduced for clar- 
ity for the two cases with central zero slope. En'or bars are the standard 
deviation of the individual measurements for each case. Each PC pixel is 
0.046". 

tion between core radii and age for a sample of LMC globular 
clusters in which the core radius seems to increase with ages 
between 1 Myr and 1 Gyr and then show s a wide range ofval- 
ues for older clusters. Using HST data, Ide Griis et alj (l2002h 
explore the matter for a sample of rich LMC globular clus- 
ters and find that young clusters tend to have small core radii 
while older clusters have an increasingly large spread of core 
radii. MAC03 explored this relationship and found that the 
relation is also valid for globular clusters around other Milky 
Way satellites besides those in the LMC. 

We concentrate in the central parts of the clusters since this 
is the region for which our technique has found differences in 
the SB shape when compared to profiles obtained from star 
counts for some clusters. Improving the measurements in this 
region and merging the results from our Galactic sample with 
those of this new sample, can help to understand their dynam- 
ical evolution. The LMC, SMC and Fornax galaxy globular 
cluster systems offer a unique window of opportunity to test 
if there are fundamental differences between systems due to 
their age. 

2. SIMULATIONS 

In paper 1 we performed a large number of simulations 
in order to establish the best method for measuring surface 
brightness profiles from HST images and also to estimate the 
uncertainties of our measurements. Results from that paper 
indicate that the only way to measure reliable surface bright- 
ness profiles from integrated light is by using high signal to 
noise images. In order to evaluate how our findings for galac- 
tic GCs translate to clusters further away, we again perform 
extensive simulations, which we describe in detail in this sec- 
tion. 
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2.1. Image Construction 

The way we create a simulated image is by adding syn- 
thetic stars on a backgr ound image using the task ADDSTAR 
in DAOPHOT ( StetsonllT98 7). The background image we use 
is a WFPC2 image of a very unpopulated field for which the 
few present stars have been cleanly subtracted. The input 
star lists are created in the same way as in paper 1. With a 
given SB profile and a luminosity function, stars are gener- 
ated randomly around a given center (the middle of the chip) 
following the two probability distributions, the surface bright- 
ness for radial distribution, and the luminosity function for 
the magnitu de distribution. The suppli ed luminosity function 
comes from I Jimenez & PadoanI (Il998h and it is corrected us- 
ing the distance modulus for the LMC. The observed lumi- 
nosity function is extended in the faint end in order to include 
unresolved background light in the simulated images, so we 
expect to recover fewer stars than the number we input. 

In paper I, we simulate SB profiles with the shape of various 
power laws. This gave us a good feel for our ability to recover 
a given central slope, but we could not test our ability to mea- 
sure turnover radius. In order to better test our method, this 
time we create a series of profiles formed by two power-laws 
joined at a break radius with a variable sharpnes s of break, 
better known as Nuker profiles dLauer et"anil995h . A Nuker 
profile is defined in the following way 



/(r) = 42- 



1 + 
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where rt, is the break radius, 4 is the surface brightness at 
the break radius, -7 is the asymptotic inner slope, -/? is the 
asymptotic outer slope, and a is the sharpness of break. By 
using these type of profiles, we are capable of reproducing 
the characteristics of observed central profiles for the sample. 
We create six different input profiles, whose parameters are 
summarized in Table [T] The radial extend of the simulated 
clusters is 400 pixels, which is equivalent to 18.4" with the 
PC pixel scale (0.046 arcsec/pixel). 

Once we have the input profiles, we proceed to create multi- 
ple realizations of a given model including different numbers 
of stars. Using various DAOPHOT tasks we add synthetic 
stars onto the background image. We use as the input point 
spread function (PSF) the one calculated for the LMC cluster 
NGC 1835 with a PSF radius of 9 pixels. Judging by the num- 
ber of recovered stars in the real data, we create images with 
three different amounts of input stars: 200,000 input stars, 
which yields ~10,000 detected stars; 50,000 input stars, giv- 
ing ^ 6,000 detected stars; and 10,000 input stars, for which 
we find ^ 2,000 stars. The majority of the real clusters in the 
sample are comparable to the first two cases. The different 
realizations have the exact same input parameters but come 
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Fig. 2. — Surface brightness profiles for three groups of simulations with 200,000 input stars. For each case (models 1, 4, and 5) the measurements from 
individual realizations (thin lines) are plotted against the input profile (thick solid line). The profiles are measured from four different images: full, subtracted, 
10% brightest stars subtracted and 3% brightest stars masked, and also from star counts. The vertical axis is on an arbitrary magnitude scale. 



from different, non-overlapping star lists. We create 10 real- 
izations for the 200,000 input stars case, and 20 for the other 
two. It is worth noting that the number of detected stars de- 
creases with increasing input central slope for the same num- 
ber of input stars. For the steepest central slope ~8,000 stars 
are found compared to the ^10,000 for the zero central slope 
cases. To avoid confusion, we always refer to the simulated 
datasets by the number of input stars rather than the number 
of detected stains. 

2.2. Center Determination 

Having an accurate estimate of the center position is a key 
step to measuring an accurate density profile. Our technique 
for finding the center of a cluster is described in detail on pa- 
per I. We take a guess center, divide the image in eight sectors 
converging at that center, count the stars in each sector and 
calculate the standard deviation of those eight numbers. We 
change to a different guess center and perform the same op- 



eration. In the end we have a grid of guess centers with a 
standard deviation value associated to them. We fit a surface 
using a spline smoothing technique (IWahballl980l : lBates et alj 
1986) and choose the minimum of this surface as our center. 

We test the accuracy of our center determination technique 
by applying it to these simulated images. Figure |2] shows 
the average measured center and the standard deviation of the 
measurements for different groups of simulations. The maxi- 
mum deviation observed is of ^1 pixels, which is equivalent 
to ^^0.3". These results are better than those in paper I. We 
believe the reason for this is that there are more stars enclosed 
in the same projected radius due to the distance difference, 
therefore the center estimation is improved. 

2.3. Surface Brightness Profiles 

We compare the results of measuring the density profile 
from integrated light versus doing it using star counts. We 
refer the reader to the detailed discussion in section 2.3 of 
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Fig. 3. — Same as previous figure but for simulations with 50,000 input stars. 

paper I about the strengths and weaknesses of each method. 
Resuhs from that paper indicate that using a robust estima- 
tor to calculate the number of counts per pixel in a given area 
is the best way to recover the central part of the profile. For 
that reason we us e the same robust estimator, the biweight 
jBeers et al.ll 19901) . for our measurements in this work. 

We measure the brightness of stars as well as their location 
from the constructed images. We use a series of DAOPHOT 
tasks to find stars and perform PSF fitting photometry with 
the same PSF that we used to construct the images. Since 
Poisson noise is included when constructing the images, this 
does not make the subtraction perfect. The process produces 
an image where all the stars have been subtracted and only the 
background light remains. We call these images 'subtracted', 
while we refer to the original unsubtracted image as 'full'. 
We measure a surface brightness profile in concentric annuli 
from the center of the cluster by calculating the biweight of 
the counts in an annulus and dividing by the number of pixels 
in that annulus. We use two sets of annuli for these measure- 



ments. The first set goes from 1 to 25 pixels in steps of 4-6 
pixels, and the second set goes from 20 to 100 in steps of 20 
pixels. The assigned radius of a given annulus is the average 
between the inner and outer radii. In paper I we note that, 
for the input profiles with steep cusps, the subtracted images 
produce a flatter central profile than the input. This is because 
crowding produces and over subtraction after the PSF fitting 
process. As done in paper I, we decide to produce alternative 
images with only the 10% brightest stars PSF subtracted and 
an image where we mask the 2% brightest stars with a radius 
of 3 pixels to try to avoid the over subtraction problem. 

Once we have the catalog of found stars we can compare 
it to the original input list and estimate how many of the in- 
put stars are found for different magnitude a radial bins. For 
the cases with 200,000 input stars We find that bright stars 
are found with an efficiency higher than 100%, meaning that 
more stars are found to be in the brightest magnitude bin than 
the number of stars that were input for that same magnitude 
bin. This effect is more pronounced in the central part of the 



SB profiles of LMC, SMC and Fornax GCs 



5 



cluster, where crowding problems are worse. The effect is 
easy to explain since one or more faint stars are likely to fall 
within the PSF disk of the bright stars and their light is mea- 
sured as if it was part of the bright stars. The difference be- 
tween input and measured magnitudes is typically ~0.1 mag, 
which is enough to push some stars from a fainter magnitude 
bin to a brighter one. For intermediate-magnitude stars, the 
same effect happens, for the regions near the core; the ef- 
ficiency for finding these stars falls to 50-70% depending 
on the shape of the profile. The efficiency for finding the 
fainter stars is lower in any radial bin; it is a few percent in 
the center and up to 50% for the regions at large radii. As 
expected, these numbers become more extreme for the cases 
with steeper central slopes, since crowding is worse then. For 
the case with 50,000 input stars, the trends are similar, but the 
numbers are less extreme. Stars in the brightest magnitude 
bin are found with an efficiency close to 100% for the cases 
with flatter central slopes. The efficiencies for the cases with 
steeper central slopes are very similar to those with 200,000 
input stars. Finally, for the cases with 10,000 input stars. The 
efficiencies for finding the input stars are all close to 100% 
except for the faintest stars in the central region of the cluster, 
which are around 70 - 80% depending on how steep the cen- 
tral slope is. The conclusion from this analysis is that when 
correction factors are calculated for star count measurements, 
the factors are dependent on the shape of the density profile 
and the number of existing stars. If one assumes the wrong 
shape or the wrong number of stars in the cluster, the correc- 
tion factors will be incorrect. 

Stars are counted in and divided by the area of each an- 
nulus. The above discussion about the efficiency for finding 
stars implies that the stars below a certain brightness are never 
found with 100% efficiency, therefore we exclude them from 
the star lists. We compare the obtained star count profiles 
with those obtained by measuring integrated light from the 
four different images (full, subtracted, 10% subtracted and 
2% masked). Results from these measurements are shown in 
Figure lZTI for the 200,000 input stars case and Figure IZTI for 
the 50,000 input case. In both figures we are showing models 
1, 4, and 5, which have central slopes of 0, -0.4 and -0.7 re- 
spectively. We find that, depending on the shape of the input 
model, the profile measured in the subtracted, partially sub- 
tracted, or masked images follow the input profile best. For 
the least concentrated cases, the measurements from the sub- 
tracted and 10% subtracted image seem to follow the profile 
best, but for the more concentrated cases, the subtracted and 
10% subtracted cases tend to look flatter in the center than 
the input profile. For these cases, the profile from the masked 
image seems to be a better choice. The star counts profiles 
are always much noisier than the light profiles in the central 
regions and they show a consistent bias in the central regions 
for the cases with steep central slopes. 

We test how well we recover the input central slope for the 
different shapes of input profiles and for the different mea- 
surement methods. Since we measure the central slope by 
taking a first derivative of the profile we need a smooth ver- 
sion of it. For this, we apply the one dimensional version of 
the spline smoother mentioned in section IZ2] This allows to 
recover information from the profile without fitting any para- 
metric model to the data. We exclude the star count profiles 
from these measurements because the central parts of the pro- 
files are too noisy for the spline smoother to get a reasonable 
fit. The first derivative of the smooth profile has a section 
toward the center where it is constant; we take this constant 



value as the measured central slope. After measuring the cen- 
tral slope for the different realizations, we calculate the aver- 
age and the standard deviation for each case. We show the 
input versus measured central slopes in Figure |23] 

We observe that the uncertainty on the slope measurements 
increases as the number of input stars decreases. For the 
10,000 input stars case, the profiles from the subtracted and 
10% subtracted images yield smaller uncertainties. These two 
cases tend to underestimate the central slopes for the concen- 
trated and rich (200,000 and 50,000 input stars) cases, while 
the slopes recovered from the masked images seem to fol- 
low the input better. For all the rich cases, the measurements 
for the model with the steepest central slope overestimate the 
slope, we think this can be due to the fact that so many stars 
are being input at the center that not enough stars are being 
input for the outer parts, which would explain the fact that we 
find fewer stars for this case. 

We test our ability to measure the input break radius by 
measuring the minimum of the second derivative, which is 
the radius at which the curvature is maximum. Our results 
show that we can measure the break radii for the simulated 
clusters to within 10% accuracy. The majority (all except two) 
of the observed clusters have a reported core radius larger than 
the one for our simulations, so we are confident that we can 
measure such break radii. 

2.4. Uncertainties 

We refer the reader to the detailed discussion in section 2.4 
of paper I about the sources for uncertainty when measuring 
surface density profiles from integrated light versus measur- 
ing it from star counts. In order to properly estimate our un- 
certainties, we compare the photometric scatter between dif- 
ferent realizations having identical input parameters with the 
biweight scatter estimate. In paper I we find that the biweight 
scatter has to be scaled in order to match the photometric scat- 
ter measured from the different realizations. For these new 
simulations we find that the scaling factors change due to the 
differences in our input shapes and of simulating clusters at 
larger distances (the number of stars on a given annulus and 
differences in PSF). As done in paper I, we compare these 
scaling factors with those obtained for r eal data from an alter- 
native method discussed in section [331 

We also estimate the error in our central slope measure- 
ments by comparing the scatter of measured slopes with the 
known input slo pe fo r every simulated cluster. The results 
are shown in Fig 12.31 We confirm what we learned from an- 
alyzing Figs 12. II and 12. II The slope uncertainties are smaller 
for the subtracted and partially subtracted cases, but they are 
biased low for the cases with steep cusps and large number 
of input stars. Also, the slope measurements are more uncer- 
tain for the clusters with 10,000 input stars. The figure sug- 
gests to take the measurements from the masked image for the 
cases with steep central profiles and the measurements from 
the subtracted or 10% subtracted for the others. In the case 
of 10,000 input stars, the subtracted case always seems to be 
better and is not biased. 

3. DATA AND ANALYSIS 
3.1. Sample 

As mentioned in paper I, there are minimum requirements 
for an image to be suitable for measurements with our tech- 
nique. The image needs to have a minimum number of counts, 
which can be obtained by having a large number of stars 
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Fig. 4. — Input versus measured surface brightness slope for the different groups of simulations. A small horizontal offset is introduced for clarity in the case of 
the two models with central zero slopes. We show the average measured slope of the individual profiles for each case. Error bars represent one standard deviation 
for all the measurements. 



present due to richness, high concentration, or by having long 
exposure times. We establish that detecting stars six mag- 
nitudes fainter than the horizontal branch with a signal to 
noise of 20 is a minimum requirement for low-concentration 
clusters. This criterion can be relaxed for highly concen- 
trated clusters and for those with a large number of stars 
{My < -7.5). Taking into account these requirements, we 
gather 30 clusters from the HST archive. The sample con- 
tains 21 clusters in the LMC, 5 in the SMC and 4 in the For- 
nax dwarf galaxy. When images are available in two filters 
(F555W and F814W), we align and combine the images in 
order to improve signal to noise. We believe we are justified 
in doing this because the color gradients for the radial range 
that we are measuring arc smaller than the photometric uncer- 
tainties. If no alternative image is available, we use the single 
F555W dataset. In general, we analyze only the chip in which 
the center of the cluster Ues, the only exception is the clus- 



ter Kron 3, for which we use all four chips. The size of one 
WFPC2 chip is large enough to contain a few core radii for ev- 
ery cluster in the sample. The scale of the CCD is 0.1 "/pixel 
for the WF chips and 0.046"/pixel for the PC chip. 

We use the WFPC2 associations from the Canadian Astron- 
omy Data Center website^ These images are spatial associa- 
tions of WFPC2 images of a given target. The raw data frames 
are processed through a standard reduction pipeline, grouped 
in associations and combined. The available data are a multi- 
group image with frames for the three WF and the PC chips. 
It is straightforward to aUgn and combine two of these im- 
ages from different filters if they belong to the same program, 
which is the case for every object with two images available 
in our sample. 

' http://cadcwww.dao.nrc.ca/ 
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TABLE 2 
LMC, SMC AND Fornax Sample. 
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3.2. Image Processing 

We process the data in the same way we do for the simu- 
lated images. We choose the frame where the cluster center is 
located, this is usually the PC frame, but for a few cases, it is 
one of the WF frames. We trim the image in order to eliminate 
the noisy edges and then proceed to perform basic photometry 
with various DAOPHOT tasks. First, we use the "FIND" task 
to make a preliminary list of detected stars, then we perform 
aperture photometry with the task "PHOT" in order to choose 
candidates for PSF construction. We find that PSF stars have 
to be chosen by hand because a single bad PSF star can have 
an important effect in the final PSF construction. Once we 
have a list of PSF stars, we perform an iterative procedure 
in which we subtract the neighbor stars to the PSF stars and 
then recalculate the PSF. In this way, the PSF construction is 
less affected by crowding. Using the final constructed PSF we 
subtract all the stars from the image, leaving behind an image 
containing only background light. We also create an image 
with the brightest 10% stars subtracted and another one with 
the brightest ^2% stars masked. 

Some of the data frames contain a small number (2 to 6) of 
very bright stars that appear saturated in the images. These 
bright stars are found at various locations on the image, but 
the closer they are to the center, the larger the effect they can 
have on our measurements by creating bumps in the profile. 
We decide to exclude these stars from our measurements by 
masking them with a larger masking radius than the one used 
for the 2% brightest stars. The clusters for which this ex- 
tra step was taken are NGC 1818, NGC 1984, NGC 2100, 
NGC 2214, and NGC 330. 



3.3. Center 

As discussed in paper I, it is crucial to measure the center 
of the cluster accurately in order to measure a reliable sur- 
face brightness profile. We use the method outlined in section 
12.21 to estimate the center coordinates for each cluster Many 
clusters in the LMC are known to have an elliptical shape. Our 
method is not affected by this as long as the cluster is symmet- 
ric in two dimensions. The first guess center is always chosen 
by eye, and the radius for counting stars is chosen so that the 
circle lies entirely inside the image, therefore, it is larger if 
the center is closer to the middle of the chip. For one cluster 
(NGC 1868), the center of the cluster is very near the edge of 
one of the chips. In this case, we calculate the center from a 
different image with lower exposure time, but with the center 
located in the middle of the chip. We report our measured cen- 
ters in table 121 We warn the reader that special care should be 
taken when using these coordinates. These centers are valid 
only using the world coordinate system (WCS) information 
contained in the header of each image. The WCS information 
for two different images can make the coordinates for same 
location vary by a few arcsec. The center is always measured 
on the primary dataset (the F555W image) when two images 
were combined. 

3.4. Surface Brightness Profiles 

Once we have measured a center, we calculate the surface 
brightness profile from the four different images of each clus- 
ter. We calculate surface brightness by estimating the bi- 
weight (as explained in section 1273] ) of the number of counts 
per unit area in a series of concentric annuli. The choice of 
the size of annuli in which we measure the profile is given 
by a trade off between spatial resolution and noise. For im- 
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Fig. 5. — Surface brightness profiles for the LMC clusters. For each cluster we show our photometric measurements (solid points), our smooth profile (solid 
line), and the EFF fit by MVM05 (dotted fine). The smooth profile comes from a fit to our photometric points inside ~ 10" and the EFF fit outside that region. 
For every panel the SB units are V mag/arcsec- . 



ages with very high signal to noise, we can use smaller steps, 
while for more sparse cases, smoother profiles are obtained 
by increasing the size of the bins at the cost of decreasing the 
spatial resolution. We use three different sets of annuli, the 
first one goes from 1-25 pixels with steps of 4-6; the second 
goes from 20-100 pixels with steps of 12-15; and the third one 
from 100-380 pixels with steps of 40-60. 

In section 123] we observe that for the simulated images the 
star counts profile tends to underestimate the profile at the 
center and it is noisier than the integrated light profile. For this 
reason, we decide not to calculate the star count profiles for 
these datasets. Also, we observe that the profiles coming from 
the unsubtracted image are always noisier than those obtained 
from the other images, so we never use the 'full' profile as 
our final result. For every set of simulations, the subtracted 
and 10% subtracted images always yield smoother profiles, 
unfortunately they show to be biased toward the center for the 



profiles with steep central slopes, so we can only use them 
when all four profiles are consistent with a central flat profile. 
If there are systematic differences between the original and 
masked profiles and the two subtracted ones in the sense that 
the first two are steeper than the latter two, then we use the 
profile from the masked image, since this is the one that traces 
the central cusp best in our simulations. 

In paper I we find the photometric zero point by integrat- 
ing our measured light profiles and comparing them to previ- 
ously obtained profiles from ground based data. We cannot 
do the same thing here because our profiles have a smaller ra- 
dial extent. For the cases in which our central profile differs 
significantly from previous measurements, the radial extent 
in which the two profiles agree is not large enough for us to 
make a meaningful comparison of enclosed light. We also ob- 
serve that the differences in shape between our measurements 
and those obtained by MAC03 are always inside the turnover 
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Fig. 5. — continued 

radius. We therefore use the data points outside the core ra- 
dius to normalize our profiles to the EFF fits by MVM05. 
We choose to normalize to these profiles because MVM05 
use the MAC03 photometric points, but they re-normalized 
them using ground based data and they correct for redden- 
ing. This brings all our measurements to a common scale on 
V,„„g/arcsec^. 

We want to make measurements of central and outer slopes, 
but our images are radially limited, so we construct radially 
extended profiles by using our measured profile inside ~ 10" 
and the MVM05 EFF fits outside that radius. We measure 
slopes by taking a first derivative of the profiles, which re- 
quires a smooth version of the profiles since noise is greatly 
amplified when taking derivatives. The smoothing is done 
by using the one dimen sional spline mentioned in section |23] 
dWahba & Wang|[l99C ft. For most profiles, there is a region 
in the center for which the first derivative is constant. We 
take this value as the central slope. For the clusters that show 



a steep central cusps, the slope sometimes changes through 
the entire radial range. In this case we take the central most 
value of the derivative as the inner slope. Since we use EFF 
fits at large radius, we expect the first derivative to reach a 
constant value outside. The measured value is expected to co- 
incide with the slope of the power-law for the EFF fits, which 
it indeed does. We take the value of the first derivative at the 
half-light radius as the outer slope. It is worth noting that the 
measured outer slope will be very different to that measured 
for Galactic clusters, since the profiles for those clusters are 
calculated from a King fit, which does not have a constant 
outer slope. 

We deproject the profiles after smoothing in order to ob- 
tain a luminosity density (LD) profile for each cluster For 
the clusters with flat central profiles we often cannot obtain 
a proper deprojection due to noise, because the noise some- 
times makes the central points be slightly fainter than the rest, 
which produces a positive slope in the smooth profile and that 
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Fig. 5. — continued 



cannot be deprojected numerically. For the cases in which 
we do obtain a proper deprojection, we measure the central 
slope of the LD distribution in the same way as we measure 
the central SB slope, by taking a first derivative. 

The traditional measurement of core radius as the radius for 
which the central luminosity value falls by half loses meaning 
when the central slope is not zero. If the central profile shows 
a cusp, then the core radius will be resolution-dependent. For 
this reason, we decide to measure what we call a break radius 
instead. The break radius is defined as the radius of maximum 
curvature, the one in which the second derivative reaches a 
minimum. This is a more systematic measure for a set of 
non-parametric profiles with different central slopes. Even 
after applying the smoothing procedure, there is still a certain 
amount of noise present in the second derivative, for this rea- 
son, we fit a high order polynomial and take the minimum of 
the fit instead of the minimum of the second derivative as our 
break radius. In paper 1 we find that the core and break ra- 
dius coincide for the clusters having a flat core, but the do not 
coincide for the cases presenting a cusp. 

The difference in shape from our measurements and the 
parametric fits will affect the measurement of the half light 
radius. Since we are using EFF fits for the outer part, and 
these fits are formally infinite, we have to truncate them in or- 
der to measure the total enclosed light. We use the tidal radius 
measured by MVM05 as a truncation radius and measure the 
half light radius for our smooth profiles. Having an estimate 
of the total luminosity and using the M/L values calculated 
by MVM05 we can estimate the total mass of each cluster 



and thus est imate the median relaxation time as described in 
iBinnev & Tremaine ( 1987) 



trh ■ 



2.06 X 10* 
ln(0.4M,/(m)) 



1/2 3/2 



We assume a mean mass of Q.5Mq as in MVM05. Results 
from these calculations are presented in table[3] 

3.5. Data Uncertainties 

We describe how we estimate uncertainties for the simula- 
tions in section 1231 The method is based on different realiza- 
tions for which shot noise from stars can be estimated directly. 
We use an alternative method to calculate the uncertainties for 
real data and we calibrate this method against that used for the 
simulations, as we did in paper 1. We assume a smooth under- 
lying stellar radial profile, so the uncertainties of the photo- 
metric points should reflect deviations from a smooth curve 
in a statistically meaningful way (i.e., have a Gaussian distri- 
bution around the mean value). We calculate the root mean 
square (RMS) difference between the smooth profile and the 
data points for the central region. The biweight yields an es- 
timate for the central location (SB value) and scale (scatter); 
this scale value is divided by the square root of the number 
of sampled pixels and used as the initial uncertainty for indi- 
vidual photometric points. We then calculate the ratio of the 
biweight to the RMS, which should represent our lack of in- 
clusion of shot noise from the stars. This ratio depends on the 
extent of the radial bins (i.e, the number of pixels used), there- 
fore we use different scaUngs for the different binnings. We 
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Fig. 6. — The same as in FigOfor the SMC clusters. 

estimate the scales for the simulations using the different real- 
izations, in order to make sure that the scalings coincide. The 
average scaling for the inner points is about 3 and about 10 for 
the outer points. These numbers are consistent with what we 
found in the simulations. Thus, we are effectively including 
shot noise from stars. The largest scalings occur for sparse 
clusters, as expected. 

In the same way as in paper 1, we calculate the uncertainties 
on slope measurements from a bootstrap technique. T he boot- 
strap approach follows that in iGebhardt et alJ (Il996h . From 
the initial smooth profile, a new profile is created by generat- 
ing random values from a Gaussian distribution with the mean 
given by the initial profile and the standard deviation from the 
photometric uncertainties. A hundred profiles are generated 
in this way and the 16-84% quartiles are measured for the 
errors. These estimated errors are compared with the scatter 
measured for the simulated cases in fig |2.3l and the two inde- 
pendent error measurements agree quite well, which gives us 
the confidence that the uncertainties calculated with the boot- 
strap method are reliable. In paper 1 we perform one more 
check on our slope uncertainties by measuring the effect of 
increasing the uncertainties on photometric points by a factor 
of two. From the bootstrap method, we find that the slope 
uncertainties increased by a modest factor, less than two, for 
most clusters. Thus, the slope uncertainties are not too sensi- 
tive to individual photometric errors. 



16 




log r (arcsec) 



4. RESULTS AND DISCUSSION 
4.1. Surface Brightness 

The measured surface brightness profiles for the entire sam- 
ple are shown in figs l3.1l to l3.1l For each cluster we show our 
normalized photometric points with error bars, and a smooth 
profile made from the combination of our photometric points 
inside ^10" and MVM05 EFF fits outside that radius. For 
comparison we show the MVM05 EFF fit and the central pho- 
tometric point used for that fit. We would like to stress that 
our measured photometric points at radii larger than ^^10" do 
not participate in the construction of the smooth fit, instead the 
EFF fits are used in that region. For about half the sample (17 
objects), the agreement between our measurements and the 
EFF fits of MVM05 is excellent, even for those cases in which 
the central photometric point by MAC03 is barely inside the 
turnover radius (such as Fornax 2) or it does not lie on top of 
the EFFfit(such as NGC 1651, NGC 1898, or2100). There is 
only one case (NGC 1754) for which our photometric points 
are fainter than the EFF fit. For the remaining 12 objects, our 
photometric points are brighter than the EFF fit by more than 
0.5 mag/arcsec^, with three objects (NGC 2019, R136, and 
Fornax 3) having differences larger than 2 mag/arcsec-. 

MAC03 identify a few clusters that they think agree with 
the expected post core-collapse (PCC) morphology by show- 
ing a power-law cusp in their central profile. NGC 2005 
and NGC 2019 are identified as clear cases of PCC morphol- 
ogy with central power-law slopes of -0.75. NGC 1835 and 
NGC 1898 are marked as good candidates for PCC morphol- 
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Fig. 7. — The same as in Fig |3.1l for the Fornax dwarf galaxy clusters. 



ogy, but they measure lower power-law slopes of -0.45 and 
-0.30 for them. Three more clusters, NGC 1754, NGC 1786, 
and NGC 1916 have incomplete profiles and are classified as 
intriguing due to their small cores, but are not placed as firm 
PCC candidates. Fornax 5 is also considered a good candi- 
date for a PCC cluster based on it's small core and central 
profile shape. Our results for these seven clusters confirm 
the presence of a steep cusp for NGC 2005, NGC 2019, and 
NGC 1916; and a shallow cusp for NGC 1786. The rest of 
the cases all show flat central cores. Our reported values of 
the central slopes are different from the power-law slopes of 
MAC03, this makes sense since they are fitting a power-law to 
photometric points on a larger radial range than that in which 
we are measuring central slopes. 

Independently of PCC morphology, we identify a few more 
clusters as having clear central cusps (with central slopes 
steeper than -0.20) such as NGC 1866, NGC 2031, Fornax 3 
and Fornax 4, and some showing shallow cusps (with cen- 
tral slopes flatter than -0.20) such as NGC 1868, NGC 2214, 
and Fornax 2. When the luminosity density central slopes are 
taken into account, a similar classification arises, NGC 1866, 
NGC 1916, NGC 2005, NGC 2019, and Fornax 3 have steep 
cusps with LD logarithmic slopes steeper than -1.00, while 
NGC 1754, NGC 1868, NGC 2031, and Fornax 4 show shal- 
low cusps with slopes between -0.2 and -1.0. The cluster 
R136 is discussed in a separate section i4.2i . We should note 
that some of these clusters have half-light relaxation times 
longer than their measured age, they cannot be expected to 
have undergone core-collapse and therefore, another mecha- 
nism has to be invoked to e x plain the central non-zero s lopes . 
iBastian & GoodwinI (l2006h : Icoodwin & Bastianl (12006 ) sug- 
gest that young star clusters can be out of Virial equilib- 
rium due to rapid gas losses, and therefore, the shape of 
their surface brightness profiles can change on relatively short 
timescales. 



4.2. R136 

R136 is known to be an extremely young object at the cen- 
ter of the 30 Doradus nebula in the LMC. It is considered to 
be a young version of a globular cluster due to its large con- 
tent of O type stars. Main sequence stars with masses a s high 
as I2OM0 have been detected in it dMassev & Hunteii ri998). 
The estimated age for the most massive stars is < 1-2 Myr 
and the mass function agrees very well with a Salpeter initial 
mass function (IMF). This makes R136 a unique and very pe- 
culiar object because it allows us to study star clusters in the 
way they looked just after formation. The surface brightness 
profile that we measure has a logarithmic central slope that is 
steeper than anything measured before for a globular cluster 
and steeper than anything predicted by dynamical models like 
core-collapse. This makes us suspect that we are not resolving 
a core or a turnover radius for this object and that our central 
slope measurement corresponds to the slope just outside the 
turnover radius for the other objects. We decide to include 
R 136 in every systematic measure we made for other clusters, 
but we caution the reader that its location in different distribu- 
tions, particularly those dealing with central SB slope, should 
be taken with a grain of salt for this reason. The central sur- 
face brightness value for this object implies a central density 
of 8 X 10^Mq/pc\ 

4.3. Combining Two Samples 

In order to explore possible correlations between physi- 
cal quantities, we combine the results for this sample with 
those for the Galactic sample from paper I. From now on, 
we will refer to the objects in the LMC, SMC, and Fornax 
dwarf galaxy as the 'satellite sample' or 'satellite clusters'. 
We compare the central slope measurements for both sam- 
ples by plotting the slope histograms side by side (Fig 14.21 ). 
We note that in both SB and LD central slopes, the satellite 
sample extends to steeper slopes than the Galactic sample. In 
total, 63% of the satellite sample is consistent with having flat 
cores, the remaining objects display a continuous distribution 
of central slopes between and -1.4 for surface brightness. 
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Fig. 8. — Histograms for surface brightness (top) and luminosity density 
(bottom) central logarithmic slopes for the LMC-HSMC-l-Fomax sample (left 
panel) and the Galactic sample (right panel). 

and between and -2.2 for luminosity density. From paper I, 
we know that 50% of the Galactic sample is consistent with 
having flat cores, a smaller fraction than for the satellite sam- 
ple. For the Galactic sample we do not find any object with 
central slopes steeper than -0.8 for SB or -1.8 for LD, we 
find two objects (R136 and Fornax 3) steeper than that in the 
satellite sample. Even when these differences are taken into 
account, the main conclusion that the slope distributions are 
inconsistent with a bimodal distribution of flat and PCC cores 
is the same for both samples. 

We plot a variety of physical quantities against each other in 
order to explore for possible correlations in both samples. The 
metallicity and age values are taken directly from MVM05. 
We observe in Figure|42]that the younger clusters, which be- 
long to the satellite sample, have a narrower metallicity and 
total mass ranges (-1 <[Fe/H]< and 3 x 10^ Mq < M,o, < 
IQ^Mq respectively) than the old ones. This can be due to 
the fact that for our sample the young clusters sample a small 
linear age regime, so they have fewer chances of populating 
the extreme mass regime. Both metallicity and total mass do 
not show any clear correlations with other physical quantities. 
The outer slope shows weak correlations for the satellite sam- 
ple in the sense that clusters with steeper outer slopes seem 
to be older, have fainter central surface brightness, and larger 
break radius. The Galactic clusters appear to have shallower 
outer slopes than the satellite ones, but this could be due to 
an effect of the difference between the Chebychev and the 
EFF fits used for each sample (see section [T4] l. The galac- 
tic clusters might show steeper outer slopes if they were ana- 
lyzed in the same way as the satellite ones, or visceversa. We 
note that there seems to be a narrow range of outer slopes be- 
tween -2 and -3 for the clusters with steep central slopes for 
both samples. There is a trend of clusters with steeper cen- 
tral slopes having brighter central surface brightness values. 
Every cluster with fi{0)< 14.0 mag/arcsec^ has a central log- 



arithmic slope steeper than -0.4. Central surface brightness 
seems to be fainter for older clusters, but this is only observed 
for the satellite sample. Regarding the break radius, we should 
clarify that the lack of clusters with break radii larger than ~4 
pc in the Galactic sample is a selection effect due to the fact 
that we required the core radius to fit on the WFPC field of 
view. Since the satellite clusters are 4 - 14 times further away 
than the average Galactic cluster, we can include clusters with 
larger break radius for the satellite sample. We note that all the 
clusters with a central surface brightness brighter than ^16 
mag/arcsec^ have break radii smaller than ~2pc. Our mea- 
sured break radius follows the sam e trend obs e rved f o r core 
radius versus age by other au thors lElson et al.l (Il989h : lElsonI 
(119921) : Ide Griis et all (|2002|) . Clusters younger tiian 1 Gyr 
have break radii smaller than 4 pc, while older clusters span 
a wide range of break radii. We notice that every cluster with 
central SB slope steeper than -0.5 has a half-light relaxation 
time shorter than 1 Gyr. Finally, the SB slope versus LD slope 
relation for the satellite clusters lies right on top of the one ob- 
served for Galactic clusters, which in turn is similar to the one 
observed for galaxies (Gebhardt et al. 1996.). 

5. SUMMARY AND DISCUSSION 

We obtain central surface brightness profiles for 21 clus- 
ters in the Large Magellanic Cloud, 5 in the Small Magellanic 
Cloud, and 4 in the Fornax dwarf galaxy. We construct and 
analyze a large number of simulated images in order find the 
most suitable way to obtain surface brightness, as well as to 
estimate our uncertainties. The profiles are constructed by 
measuring integrated light with a robust statistical estimator. 
We combine HST/WFPC2 images in two filters (F555W and 
F8 14W) when available and present profiles normalized to V- 
band magnitudes. 

When our results are compared with previous results that 
use different analysis techniques, we find very good agree- 
ment for ^60% of the sample. For the remaining 40%, our 
central photometric points are brighter than previous measure- 
ments. Most central surface brightness values change from 
previously reported ones with values up to two magnitudes 
brighter. For some objects in the sample, the new measured 
surface brightness profile is no longer compatible with a flat 
core parametric fit. The main reason for this difference is the 
increased spatial resolution of HST, but also because we use 
a non-parametric estimate as opposed to the traditional King 
model fits. For some of the observed profiles the departures 
from a flat core model are small, but significant. We confirm 
the existence of a steep central cusp for three clusters previ- 
ously classified as post core-collapse. We also find a subpop- 
ulation of objects with shallow cusps with logarithmic central 
slopes between -0.2 and -0.5. When we plot a variety of 
physical quantities searching for correlations, we find indi- 
cations that the younger clusters tend to have smaller break 
radius, shallower outer slopes, and brighter central surface 
brightness. In particular, the youngest cluster in the sample, 
R136, shows the steepest central profile and the brightest cen- 
tral surface brightness. We also observe a clear correlation in 
which the clusters with the steepest central slopes are the ones 
with the brightest central surface brightness. 

There have been two mechanisms explored for producing 
cusps in star clusters: core-collapse and the presence of an in- 
termediate mass black hole in the center of the cluster. A de- 
tailed discussion and references on this subject can be found 
on section 1.2 of paper 1. The range of 3-dimensional density 
slopes is wider for core-collapse than for black hole models. 
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Fig. 9. — Surface brightness central and outer logarithmic slopes, logarithmic break radius (in parsecs), central surface brightness, metallicity, logarithmic 
age, and total mass plotted against each other for the LMC+SMC+Fornax sample (open points) and the Galactic sample (solid points). R136 is shown as a star 
symbol. We also show on the top right comer two panels with SB slope versus half-light relaxation time and versus LD slope (the solid line represents 'LD slope 
= SB slope + 1'). The distances to the clusters are assumed to be 45 kpc for the LMC, 60 kpc for the SMC and 140 kpc for the Fornax dwaif galaxy. 



but they both center around the same value, ~ -1.65. How- 
ever, only the four clusters with the steepest profiles in our 
sample fall in this range. In the case of core-collapse, the 
slope depends on the mass of the stars used to measure the 
profile and of those that dominate the mass of the core, so this 
could extend the range toward shallower slopes. Another fac- 
tor of uncertainty is the time dependence of the core-collapse 
model when the core goes through gravothermal oscillations. 
According to Fokker-Planck simulations, a star cluster will 
spend a considerable amount of time in between successive 
collapses, where the light profile resembles a King model with 
a flat core. Unfortunately, these models do not give enough 
details about the slope of the density profile during intermedi- 
ate stages of post-collapse bounce, or about the time spent on 
intermediate stages. 
As discussed in paper I, an alternative way to explain the 



shallow cusps, is by invoking the result by Baumgard t et alj 
(2005). They perform simulations of stellar clusters with in- 
termediate mass black holes in their centers and find that, after 
a Hubble time, the projected density distribution of the clus- 
ters show shallow cusps with slopes around -0.25. For this 
sample, we find 4 objects that fall within this regime, but only 
kinematical data can confirm the possible existence of a cen- 
tral black hole for these objects. 

The observed co rrela tions with age observed for the Satel- 
lite sample (section l431 l point out to the possibility of clusters 
having very concentrated profiles during early stages of their 
evolution. In particular the break radius-age relation observed 
here and by many authors tells us that the size of cores de- 
pends on the dynamical evolution of clusters. The input den- 
sity profiles for various dynamical simulations have almost 
always been characterized by King or Plummer models. This 



SB profiles of LMC, SMC and Fornax GCs 



15 



TABLE 3 
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0.14 
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could be biased toward large flat cores, when more concen- 
trated profiles could be more appropriated. This is true for 
core-collapse models, as well as models containing a central 
black hole. 

Tables with the complete photometric points and smooth 
profiles for every object in this sample can be found in the 
CDS- Vizier Service. 
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